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Entanglement sharing among pairs of spins in Heisenberg antiferromagnets is investigated using 
the concurrence measure. For a nondegenerate S=0 ground state, a simple formula relates the 
concurrence to the diagonal correlation function. The concurrence length is seen to be extremely 
short. Various finite clusters are studied numerically, to see the trend in two dimensions. It is 
argued that the concurrences vanish between pairs of spins which are not nearest neighbors, for 
the linear chain and square lattice antiferromagnet ground states. For the triangular lattice and 
Kagome' lattice, nearest-neighbor concurrences also vanish. The concurrences in the maximal-spin 
states are explicitly calculated, where the concurrence averaged over all pairs is larger than the S—O 
states. 



I. INTRODUCTION 

In recent years quantum entanglement has emerged 
as a common platform for scientists working in various 
fields such computer science, physics, mathematics and 
chemistry In particular, quantum entanglement of 
spin-^ degrees of freedom, qubits, has been studied ex- 
tensively, due to their importance for quantum comput- 
ers, not to mention their well-known applicability in vari- 
ous condensed-matter systems, optics and other branches 
of physics. For a pure state of many qubits, quantum 
entanglement, which is quantified by the von Neumann 
entropy of the reduced density matrix, is a measure of 
how a subsystem is correlated to the rest of the system. 
The key ingredient for entanglement is the superposition 
of states. A linear combination of two pure states, both 
of which are entangled, viz. having a non-zero entropy 
for a subsystem, can exhibit no entropy at all. And the 
reverse situation is also possible, two pure states with no 
entanglement, can give rise to entanglement on superpo- 
sition. 

For a system consisting of a large number of qubits, 
how different pairs share entanglement in a pure state 
cannot easily be specified, even with known diagonal and 
off-diagonal correlation functions. We will show below, 
for the ground state Heisenberg antiferromagnet with 
zero total spin, the concurrence measure can be speci- 
fied in terms of the diagonal correlation function alone. 
In a pure state of many qubits, a subsystem of two qubits, 
in general, will be in a mixed state. A mixed state den- 
sity matrix can be written as a decomposition over pure 
state density matrices, with a large number of possible 
decompositions over pure states. The entanglement of 
the pair of qubits is the average entanglement of a de- 
composition, minimized over all possible decompositions. 
Starting from a given many-qubit state, with a density 
matrix p = \^ >< -01, the reduced density matrix Rij 
for a pair of qubits is constructed by performing a par- 
tial trace over the rest of the qubits to be eliminated, 
Rij = trp. In general the reduced density matrix rep- 
resents a mixed state for the pair of sites labeled 
The von Neumann entropy calculated from the eigenval- 
ues r„ of Rij , as — ^ r„ log2 r„ quantifies entanglement 



of this pair with the rest of the qubits. The concurrence 
measure 2] has the important information as to how these 
qubits are entangled among themselves, and is given as 

^ \l/2 \l/2 \l/2 \l/2\ n\ 

Cij — max(0, Aj' — \^ - X^' — A4' ). (1) 

In the above are the eigenvalues in decreasing order 
of the matrix where R is the time-reversed matrix, 

R — (7y X (TyR^dy X (J y . 

The two-site concurrence depends only on the struc- 
ture of entanglement of a many-qubit state without ref- 
erence to a Hamiltonian. However, we would like to in- 
vestigate the pairwise entanglement in the ground state 
of the Heisenberg Hamiltonian. We will study the entan- 
glement in the ground state Heisenberg antiferromagnet, 
with a Hamiltonian 

^ = E ^ ^^tSj + Si S+) + 5f 5| , (2) 

The sum is over all nearest-neighbor pairs on a given lat- 
tice. Working in a diagonal basis of for every site i, 
there are two states per site, viz., t, and [. A many- 
spin state can be characterized by the number of down 
spins, as the total z-component of the spin, 5^ = ^ Sf , 
is a good quantum number, in addition to the total spin 
quantum number S. The above Hamiltonian is also in- 
variant under time reversal. The ground state will belong 
to S'^ = subspace, as this sector has a representation 
state from every spin sector labeled by S — 0, 1, ..N/2. 
For most lattices, and finite clusters the ground state be- 
longs to S' = 0, for example the linear chain, the square 
lattice, the triangular lattice, Kagome' lattice, and fi- 
nite clusters that will be considered below. However, the 
ground state energy and the wave function are known 
only for the case of the linear chain through the Bethe 
Anzatz The highest-energy state has S = N/2, which 
is the maximum value for the total spin, with a degen- 
eracy of iV + 1 corresponding to the different values of 
~ ~N/2, ..N/2. This maximal-spin state would be 
the ground state, if the interaction is ferromagnetic in 
the above Hamiltonian (by changing the overall sign of 
the Hamiltonian). 
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II. 



CONCURRENCES IN 5 = STATES 



Let us consider a state with m down spins 

\lp >=^ll){Xi..Xm,)\xi..Xm > (3) 

where the ordered labels xi..Xrn denote the locations of 
the down spins. Let us rewrite the state as 



\tp>= ^ iP{si..sn)\si..sn > 



(4) 



where Si is the eigenvalue of 5? for the i'th qubit. The 
reduced density matrix Rij of two sites i and j {j > i) 
has the matrix elements 

-Rsiisj = X] '4'*{si,--Si-s'j..Sn)tp{si..Si..Sj..Sn) (5) 



where the sum is over all spin variables except at sites 
i and j. Since the many-qubit state has a definite eigen- 
value for 5^, we have [Rij. Sf + 5|] = 0. This in turn 
implies the following structure for the reduced density 
matrix 



Ri 



* \ 



vij 

wiij z*j 

Zij W2ij 

Uij j 



(6) 



In the above we used the two-qubit basis states | \\> 
,1 Ti>,l iT>,l ii>, where | t>,| i> stand for eigen- 
states of SI with eigenvalues 1/2,-1/2 respectively. 

The matrix elements can be expressed in terms of the 
expectation values, < A >=< tlj\A\il) >, as 



1 



,1 



1 



,1 



(7) 



Zi,j =< S+S- >, 



(8) 



(9) 

The diagonal matrix elements are simply related to the 
diagonal correlation functions, and the off-diagonal ma- 
trix element is just the off-diagonal correlation function. 
The concurrence for the two sites has now a simpler form 



Cij = 2 max (0, \zij \ - s/uijVij)- 



(10) 



As can easily be seen from above, the concurrence mea- 
sure uses both diagonal and off-diagonal correlation func- 
tions. Whether or not two sites have a non-zero con- 
currence is not at all intuitive, given a specific state 
with known correlation functions. For a long-ranged 



concurrence, the necessary condition is the existence of 
off-diagonal long range order (ODLRO). If there is no 



off-diagonal long range order (ODLRO), 



0, as 



l^i — TjI 00, the concurrence of two sites far apart 
would go to zero. Thus, the existence of ODLRO is a 
necessary condition for a long-ranged concurrence. How- 
ever, even with a ODLRO, long-ranged concurrence may 
be absent, if |zy | < ^UijVij. We will see below that the 
pairwise concurrence is short ranged in the ground state 
of Heisenberg spin systems. For a system with a large 
number of qubits, only the nearest-neighbor concurrence 
seems to be nonzero. 

Let us consider a nondcgencratc 5 = state, which can 
be represented in — subspace. Because of the time- 
reversal symmetry of the Hamiltonian, the wave function 
can be chosen to be real. This would imply that the 
off-diagonal matrix element Zij of the reduced density 
matrix will be real. Further, the rotational symmetry 
of the Hamiltonian would imply, for the S = sector, 
< SfS^ >=<SfS] >=< SfS^ >. Exploiting this prop- 
erty and denoting the diagonal correlation function by 



^ii =< Si 



5f >, we have 



< StS7 >= 2r, 



(11) 



The time-reversal invariance also implies, < >= 0, 
in this state. This will simplify the diagonal matrix ele- 
ments of Rij . Thus the diagonal elements of the reduced 
density matrix are 



(12) 



In addition the diagonal correlation function has the 
property \rij\ < 1/4, since we are dealing with spin- 1/2 
species. Now the concurrence can be specified completely 
in terms of the diagonal correlation function alone. 

The diagonal correlation function between two arbi- 
trary spins can take both positive or negative values. 
For a pair of qubits with Tij > 0, Cij = 2max(0, 2rij — 
J — Tij) = 0. The above argument would imply that the 
concurrence is zero for a pair of qubits, if their diago- 
nal correlation function is positive, in a nondegenerate 
5 = state. On a bipartite lattice, if two spins belong to 
the same sublattice. them the diagonal correlation func- 
tion would be positive, implying the concurrence here is 
zero. On the other hand, when Tij < 0, then we have 
y^UijVij = 1/4— \Tij\, so that the concurrence is given by 
Cij = 2max(0,3|rij| — i) = 6(|ry| — j^) if positive, and 
otherwise zero. This gives us a simple test for a nonzero 
conciuTcnce, viz. if Tij < 0, iFjjl > 1/12. Thus in a 
nondegenerate 5 = state, the concurrence for a pair 
of spins is specified entirely by their diagonal correlation 
function, as 



Cij = 0, for Tij > 

= 0, ioiTij <0,\Tij\< 



12 



= 6(|r,,|-l), forr,, <o,|r,,|>l (13) 
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FIG. 1: A Tic-Tac-Toe cluster of twelve atoms. There are 
eight inequivalent pairs here. Only two inequivalent pairs of 
sites have nonzero concurrences in the ground state: Ci,2 = 
0.08, Ci,4 = 0.41. 



The above formula substantially simplifies the procedure 
of calculating the pairwise concurrences in a many-qubit 
S — state. We simply calculate the diagonal spin-spin 
correlation function, and the concurrence is completely 
specified from the formula outlined. In the following sec- 
tion, we will effectively use the formula to study the en- 
tanglement distribution and sharing among various pairs 
in finite clusters and higher dimensional lattices. 

Let us turn now to computing the actual values of 
the concurrences in one-dimensional systems. Now, for 
a closed chain of 12 sites, the above implies, denot- 
ing Cij = Cr where r = j ^ i, in the ground state 
Ci = 0.398, and = for r > 1, where we have used 
the numerically computed diagonal correlation functions, 
Ti = -0.1496, = 0.0626, = -0.0553 etc.. As the 
size of the chain is increased, the magnitude of the corre- 
lation function Tr will take the limiting value from above. 
This would imply Cr = 0, for r > 1 for an infinite chain. 
Ci can be estimated from the ground state energy in the 
thermodynamic limit. The ground state energy of the 
Heisenberg Hamiltonian in general can be written as 



E„ 



e„N = 3iV„ri 



(14) 



where eg is the ground state energy per site, and Nn 
the number of nearest-neighbor pairs of spins. Thus the 
nearest-neighbor concurrence can be written as 



N 



(15) 



when positive, and otherwise zero. For a linear-chain 
antiferromagnet ground state, a similar formula has been 
derivedp|. The ground state energy is known from Bethe 
Ansatz, = In 2 — 1/4, and — N, which determines 
the nearest-neighbor concurrence, Ci — 2 In 2—1 « 0.386. 



III. FINITE CLUSTERS AND EXTENDED 
TWO-DIMENSIONAL LATTICES 



In two-dimensional systems there are no analytical re- 
sults for the ground state of Heisenberg antiferromag- 
net systems, like the Bethe-ansatz solution for the linear 
chain. Here, we will take recourse to numerical diagonal- 
izations of finite clusters, and use variational estimates 
for the ground state energies on various extended lat- 
tices. The main difference between two-dimensional sys- 
tems and the linear chain considered above is the num- 
ber of nearest neighbors, apart from various local lattice 
structural differences. To gain insight into the entangle- 
ment distribution in two dimensions, we will numerically 
study below various finite clusters, which are drawn from 
square, triangular and Kagome' lattices. 

Let us first consider a toy problem of a Tic-Tac- 
Toe cluster of twelve sites, as shown in Fig.l. This 
cluster is drawn from a square lattice. Here there 
are two inquivalent sites, site 1 is equivalent to the 
other sites on the boundary, and site 4 is equivalent 
to the other sites forming the vertices of the inside 
square. There are eight inequivalent pairs of sites: (1,2), 
(1,3),(1,4),(1,5),(1,9),(1,10),(4,5), (4,9). From the nu- 
merical exact diagonalization for the S — ground state, 
only two inequivalent pairs have nonzero concurrence: 
Cu2 = 0.08, CiA = 0.41, and Ci,3 = Ci,5 = Ci,9, Ci,io = 
^4.5 — 6*4,9 = 0. In this finite cluster, the third-neighbor 
bonds, (1,2), (1,6), (1,7), (1,11) have concurrence equal to 
0.08, where as the second-neighbor bonds (1,5), (1,8), (1,3) 
etc. have zero concurrence. This can be ascribed to 
finite-size effects, and to the fact that this cluster does 
not have the same symmetry as the square lattice. 

Let us now consider a 16-site square lattice cluster, 
with periodic boundary conditions. Here all sites are 
equivalent, with coordination number equal to four as on 
an infinite square lattice. This cluster has also the bipar- 
tite structure of the square lattice. There are only four 
inequivalent pairs, due to the periodic boundary condi- 
tions. For a given site, the second and third neighbors 
sit on the same sublattice, implying that the diagonal 
correlation function is positive, F — 0.071 from the nu- 
merical calculation, and hence a zero concurrence. The 
numerical values of the correlation function are -0.117 
and -0.067 for the nearest and furthest neighbor pairs 
respectively. Only the nearest-neighbor concurrence is 
nonzero, Ci — 0.202 which is almost half the value we 
got for the linear chain. To study larger clusters, one has 
to increase the number of sites to 24, to keep the bipartite 
structure intact. Then exact numerical diagonalization is 
quite difficult for such large clusters. 

For the infinite square lattice the exact ground state 
energy is not known from analytical solutions, however, 
there are excellent estimates of upper bounds on Cg. 
For the square-lattice ground stated, we have an up- 
per bound estimate \eg\ ~ 0.66. This in turn, along with 
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FIG. 2: (a) A triangular-lattice cluster with ten atoms. Here all the sites are not equivalent. There are four inequivalent sites 
labeled 1, 2, 4, 5. (b) A 16-site triangular cluster, with periodic boundary conditions. The extra nearest-neighbor bonds arising 
from the boundary conditions are shown with dotted lines. Here all sites are equivalent, with six nearest neighbors, like in the 
infinite triangular lattice. However, the horizontal bonds (sixteen in number) are not equivalent to the slanted bonds (32 in 
number). 



Nn = 2N, implies 

Ci = ~ ^ ~ 0-16 square lattice. (16) 

And the next-neighbor concurrence is zero, as the two 
spins sit on one sublattice of the square lattice. We 
conjecture that Cr = 0,r > 1 for square lattice also. 
This is because the concurrence reduces when the num- 
ber of neighbors increases, as is evident from the nearest- 
neighbor concurrence (it reduced from 0.386 to about 
0.16, as the number of nearest neighbors increased from 
2 to 4, from the linear chain to the square lattice). And 
in conjunction, we have the fourth-neighbor concurrence 
zero from the sixteen-site cluster studied above. 

Let us now turn our attention to the triangular lat- 
tice, which has an increased number of nearest neigh- 
bors, along with a new ingredient of frustration among 
the spins on a triangular block, for minimization of en- 
ergy. Again we take recourse to numerical diagonaliza- 
tion of finite clusters, in the absence of exact results for 
the ground state Heisenberg antiferromagnet. We will 
consider two different clusters here. Let us first look at 
a ten-site cluster shown in Fig. 2a, which is drawn from 
a triangular lattice. Though all sites are not equiva- 
lent here, but it has the basic triangular and hexago- 
nal block structure of a triangular lattice. From the nu- 
merical exact diagonalization for the ground state, only 
three inequivalent pairs have nonzero concurrence. The 
horizontal nearest-neighbor bond on the boundary has 
a concurrence, Ci,2 = 0.26, along with the three other 
equivalent bonds. The slanted nearest- neighbor bonds 
on the boundary have a concurrence, Ci^4 = 0.347. The 



bond at the center has a concurrence, C^^ = 0.336. All 
the other nearest-neighbor bonds have zero concurrence, 
along with all other pair concurrences between sites fur- 
ther apart. The magnitude of the nearest-neighbor pair 
concurrence has reduced here from the linear chain on 
account of the increased number of neighbors and the 
frustration associated with the triangular blocks. 

Let us now consider a 16-site triangular cluster with 

periodic boundary conditions, shown in Fig. 2b. The sites 
on the boundary get extra neighbors, due to the bound- 
ary conditions, indicated by dotted lines in the figure. 
For instance, site labeled as 1 has extra neighbors in sites 
labeled as 4,13 and 14, and so on for the other sites on 
the boundary (see figure). Now all the sites are equiv- 
alent, with six nearest neighbors as on an extended tri- 
angular lattice. This cluster has the ingredients of an 
infinite triangular lattice: six coordination, triangular 
block structure and associated frustration, equivalence 
of all the sites. One important ingredient that is missing 
here is the tripartite structure. There are two inequiva- 
lent nearest-neighbor bonds here, horizontal and slanted 
bonds. Starting from any site and traversing four hori- 
zontal bonds, one can return to the original site, whereas 
by traversing four slanted bonds one cannot return, as 
can be seen from the figure. From the numerical diag- 
onalization, we get for the concurrence for the 16 hor- 
izontal bonds, Ci^2 = 0, and for the 32 slanted bonds, 
Ci,5 = 0.05. The concurrence is zero for all other pairs 
of sites further apart. Here, the concurrence is furthur 
reduced from that of the square-lattice clusters due to 
increased neighbors and frustration. This discrepancy 
between the horizontal bonds and the slanted bonds will 
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disappear as the number of sites is increased, and the 
tripartite structure restored for an extended triangular 
lattice. 

For an extended triangular lattice, we will use the vari- 
ational estimates for the ground state energy, for calcu- 
lation of nearest-neighbor concurrences. For a triangular 
lattice, we have Nn = 3iV, which implies 

le I 1 

Ci = - -) for triangular lattice (17) 

if positive, and otherwise zero. The upper bound esti- 
mate for the ground state energy|^ Q is \eg\ « 0.53. 
Hence, the nearest-neighbor concurrence for the triangu- 
lar lattice is zero. The next-neighbor concurrence cannot 
be argued to be zero based on a bipartite structure, as 
we did above for the linear chain and the square lat- 
tice. However, we can argue that the absolute value of 
the correlation function would decrease as the separa- 
tion increase. Also, from the 16-site cluster with six co- 
ordination discussed above, other than nearest-neighbor 
concurrences are vanishing. We expect as the size is in- 
creased, they will still be zero. This would imply Cr = 
for all separations on a triangular lattice. 

We would like to study the Kagome' lattice clusters 
now. The obvious cluster that can be considered is a 
David star cluster of twelve sites shown in Fig. 3, drawn 
from a Kagome' net. This cluster too has triangular 
blocks, along with the associated frustration, apart from 
unfilled hexagonal blocks unlike the triangular lattice. 
Numerically it is seen that this cluster has two degener- 
ate S — ground states. The pair concurrences in one of 
the ground states, for the nearest-neighbor bonds can be 
calculated from the reduced density matrices. The only 
nonzero nearest-neighbor concurreces are for alternate 
bonds on the boundary, traversing clockwise, we have 
Ci4 = C5^7 = Cii,io = ^^12, 9 = Cs,6 = ^2,3 = 0.55. Sim- 
ilarly the other ground state has nonzero concurrences 
only for pairs (1,3), (2,6), (8,9), (12,10), (11,7), (5,4). 
This is because these ground states have a large ampli- 
tude for singlets for these pairs. This can be seen eas- 
ily by viewing the cluster shown in figure as a closed 
chain with every alternate site having a next-neighbor 
interaction. If every site has a next-neighbor interaction 
whose strength is half that of the nearest-neighbor in- 
teraction, then this reduces to Majumdar- Ghosh model, 
which has a two-fold degenerate valence-bond ground 
state j7|. But here, the next-neighbor interaction strength 
has been transfered to every second next-neighbor bond. 
But, still, the amplitude of the valence-bond configura- 
tion is predominant, giving rise to dominant contribution 
for concurrences only for alternate bonds. All the other 
pairs have zero concurrence. However, since the ground 
state is degenerate, one can change these concurrences 
by suitable linear combinations of the degenerate states. 
Larger clusters drawn from the Kagome' lattice, keeping 
the ingredients intact, are too large for exact diagonal- 
ization. Again, we will turn to the variational estimates 
for the ground state energy, to infer the nearest-neighbor 
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FIG. 3: A Kagome'-lattice David star cluster with twelve 
atoms. The ground state here is two-fold degenerate, with 
large amplitude for valance-bond configurations. 



concurrences. 

For the case of an extended Kagome' lattice, we 
have Nn = 27V, and the estimate for the ground state 
energy0,|3|, \eg\ « 5/12, would imply that the nearest- 
neighbor concurrence is zero. Analogous to the triangular 
lattice, we can argue Cr = for all separations. Though 
the square lattice and the Kagome' lattice have the same 
number of nearest neighbors, Ci is zero for Kagome lat- 
tice. This we can attribute to the frustration present 
in the Kagome' case, due to the presence of triangular 
blocks 0. 

We can increase the number of nearest neighbors by 
considering a c?-dimensional hypercubic lattice, with 2c? 
nearest neighbors for every site. For large d, the ground 
state energy has an upper bound \eg\ w d/4, and Nn = 
Nd which will yield C = for the nearest neighbors. 
Similarly we expect all pairs of sites to have zero con- 
currence. This is consistent with our argument, that an 
increase in the number of neighbors will result in a re- 
duction in the nearest-neighbor concurrences. We can 
also increase the effective number of neighbors by includ- 
ing a long-ranged interaction, by making the sum in the 
Heisenberg Hamiltonian (Eq.2) run over all pairs. Then 
each qubit interacts with all other qubits by the same 
strength. This gives usH ^Y.Si- Sj^ {S^ - 37V/4)/2, 
the eigenvalue depending only on the total spin quantum 
number S. A larger value for the total spin will lead 
to a larger energy. The lowest eigenvalue is for 5 = 0, 
eg = -3/8. Now, we also have CgN = 3N{N - l)T,j/2, 
as all the two-point correlators would be identical, and 
Si ■ Sj = 3T,j for a S" = state. This yields C = 0, for 
all pairs of qubits. 
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IV. MAXIMAL-SPIN STATES 



the function are C{S = N/2, S'^ = ±N/2) = 0, and 



Let us turn our attention to the maximal-spin states, 
with S — N/2, which are the ground states for a fer- 
romagnetic interaction, and are maximal-energy states 
for the antiferromagnet, with an energy Emax = N/A. 
The structure of the maximal-spin states will be same 
for all clusters and lattices, irrespective of coordination 
number and local structure. There are TV + 1 such 
states, with one state each in a subspace with a given 
value 5*^ = -N/2,. .N/2. Let us write 5^ = N/2-m 
where m is the number of down spins in each basis state, 
m = 0, 1..-/V for various subspaces. It is straightforward 
to write down the eigenstate for any lattice as 



\S 



N 



N 



m >= 



1 



\^> 



(18) 



where the sum is over basis states each with m 

down spins exactly. Since every basis state is given 
the same weight, the reduced density matrix for any 
pair of spins will have the same matrix elements, and 
the correlation functions, both the diagonal and the 
off-diagonal, will be same for every pair. This implies 
Cij = C for all the Np = N{N - l)/2 pairs. The diag- 
onal correlation function, thus, can be written as Tij = 
i^Np) Eaiipairs < >= (^^' " N/4)/N{N-l). For 

values of ~ 0{N), there is a long-ranged order. This 
would cause a decrease in the concurrence for a given 
pair, with a concurrence ~ 0{1/N). However, the av- 
erage concurrence for this case is more, as we will see 
below, than the concurrence in the S = state, that we 
discussed for various lattice systems above. Thus, the 
off-diagonal matrix elements of Rij are evaluated as 



N N 2 



(22) 



The concurrence is maximum when the number of 
down spins is either 1 or N-1. Though all pairs have 
the same concurrence, however, the concurrence is not 
long ranged, since in the thermodynamic limit, N oo, 
the concurrence goes to zero as 1/iV. But, in comparison 
the state with S — N/2, — has a larger average con- 
currence than the state with S = = that we studied 
earlier, where the best average concurrence for the lin- 
ear chain is < C >w 0.8/(A^ — 1). Thus, increasing the 
spin from S = 0, for the ground state, to 5* = N/2, 
for the maximal-energy state of the antiferromagnet, the 
nearest-neighbor concurrence has decreased from 0(1) to 
0{1/N). At the same time, the average over all pairs has 
improved, implying a better entanglement sharing in the 
maximal-spin states. Thus, the maximal-spin states with 
S = N/2, \ S^\ < N/2 have a better entanglement shar- 
ing among the qubits, with the best sharing in the states 
S = N/2, \ S^\ = N/2 — 1, with an average concurrence 
< C >= 2/N. The sectors with S = N/2 - 1, N/2 - 2, 
namely the one-inagnon and two-magnon states, have 
been investigatedflOl However, it remains to see if 
states with intermediate values of energy and/or with in- 
termediate values of spin < 5 < N/2 — 2, can exhibit 
more average concurrence. Most important states are the 
first excited sates of the antiferromagnet, with 5 = 1, 
where the concurrence results are yet to be known. It 
would be interesting to investigate the concurrence as a 
function of S, in all spin sectors. 



(t-^^)( 



1) 



(f + 



N{N- 1) 



N{N- 1) 



(19) 



Now, for the off-diagonal matrix element Zij =< 
Sj'S^ >, the contributions come from all configurations 
with I at site j and t at site i. There are N — 2Cm-2 
such configurations. Hence, 



(f 



N{N- 1) 



(20) 



Putting all things together, the concurrence for any pair 
of qubits is, with S = N/2, 5^ = N/2 - m. 



N{N- 1) 



m,[N — to) 



(21) 

A similar formula has been given in 9] . Some values of 



V. CONCLUSION 

We have investigated the quantum entanglement shar- 
ing in Heisenberg antiferromagnets using the concurrence 
measure. For S=0 unique ground states, a simple formula 
relates the concurrence between an arbitrary pair of spins 
to their diagonal correlation function. This substantially 
simplifies the calculation of concurrences, both numeri- 
cally and analytically. The nearest-neighbor concurrence 
is directly related to the ground state energy. For larger 
number of nearest-neighbors, the concurrence is smaller, 
as has been seen from studying a number of finite clus- 
ters, and from the study of large-d hypercubic lattice and 
longer-ranged interactions. It has been argued that other 
than nearest-neighbor pairs have zero concurrence in the 
ground state for the linear chain and the square lattice. 
It is shown that the nearest-neighbor concurrences also 
vanish in the ground state of triangular and Kagome lat- 
tice antiferromagnets. For maximal-spin states, explicit 
formulas are given for the pairwise concurrences, and the 
states with S — N/2, = N/2 — 1 show maximum en- 
tanglement sharing. 
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